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Spin-dependent electronic transport through multiferroic Co/PbTiOs/Co tunnel junctions is stud¬ 
ied theoretically. Conductances calculated within the Landauer-Biittiker formalism yield both a 
large tunnel magnetoresistance (TMR) and a large tunnel electroresistance (TER). On top of this, 
we establish a four-conductance state. The conductances depend crucially on the details of the 
electronic structure at the interfaces. In particular, the spin polarization of the tunneling elec¬ 
tronic states is affected by the hybridization of orbitals and the associated charge transfer at both 
interfaces. Digital doping of the PbTiOs barrier with Zr impurities at the Ti02/Co2 interface 
significantly enhances the TMR. In addition, it removes the metalization of the barrier. 

PACS numbers: 72.25.Mk, 73.40.Rw, 75.85.+t 


I. INTRODUCTION 

Recently, the functionality of magnetic tunnel junc¬ 
tions has been enhanced by replacing the insulating bar¬ 
rier, typically a band insulator such as MgO (e. g. Ref.0), 
by a ferroelectric. This combination of ferromagnetism 
in the electrodes and ferroelectricity in the barrier leads 
to a multiferroic tunnel junction (MFTJ) in which tun¬ 
nel magnetoresistance (TMR) and tunnel electroresis¬ 
tance (TER) show up simultaneously. By independently 
switching the mutual orientation of the magnetizations in 
the electrodes—say, from parallel (P, to antiparallel 
(AP, ff )—and reversing the ferroelectric polarization in 
the barrier—say, from left («—) to right (—>>)—the tunnel 
conductance may take four different values, leading to a 
nonvolatile four-state memory device^. 

A four-conductance state has been reported recently 
for Co/PbTiOs/LaSrMnOs tunnel junctions^, in which 
a perovskite barrier of 3.2 nm PbZro.2Tio.8O3 (PZT) has 
been epitaxially grown on the ferromagnetic and almost 
half-metallic perovskite La(Sr)MnO3(001) (LSMO). An 
intriguing observation is a sign change of the TMR: 
switching the ferroelectric polarization in PZT from 
pointing towards LSMO to pointing towards Co reverses 
the TMR from +4 % to -3 %. 

For magnetic tunnel junctions, e. g. Fe/MgO/Fe, it has 
been shown that details of the interfaces play a decisive 
role for the TMR£~— : their geometry and their electronic 
structure govern the spin-dependent transport^. Con¬ 
cerning Co/PZT interfaces, ab initio calculations^ sug¬ 
gest that (i) the interface is Ti(Zr)02-terminated. More 
precisely, the interfacial Co atoms are placed in line with 
O. (ii) Reversal of the ferroelectric polarization of PZT 
may also reverse the spin polarization at the interface. 
It is conceivable that this switching is responsible for the 
experimentally observed sign change of the TMR. The 
interface between LSMO and PZT could be Pb0/Mn02 
since the ferroelectric films grow epitaxially and in com¬ 
plete unit cells on clean substrate s 10 ! 11 ; this scenario is 


supported by first-principles calculations^. 

To explain the origin of the inverted TMR4£ in 
Co/PZT/LSMO, Pantel et al. point out the role of 
LSMO (Ref. 0)- A spin polarization of the LSMO 
surface as large as 95 % has been deduced from 
transport measurements^. Thus, both interfaces of 
Co/PZT/LSMO may contribute significantly to the 
above-mentioned peculiarities of the TMR. In Refer¬ 
ence 0 a contribution from resonant tunnelingi^ - — via 
electronic states localized within the barrier has been 
ruled out. We recall that resonant states can alter the 
spin polarization of the tunneling electrons^. 

In this work, we investigate spin-dependent tunnel¬ 
ing in Co/PbTi0 3 /Co and Co/PbTi0 3 -Ti(Zr)0 2 /Co 
MFTJs by means of first-principles electronic structure 
and transport calculation s 17 ! 18 . These tunnel junctions, 
with their two Co electrodes, have been chosen to exclude 
effects of an LSMO electrode. We focus on the electronic 
states at the Co/PTO interfaces (PTO = PbTi0 3 ) and 
address in particular how the interfacial magnetoelectric 
coupling affects the local magnetic moments. For this 
purpose, we compare results for tunnel junctions with 
one Co/Ti02 or one Co/Zr02 interface while keeping a 
Pb02/Co interface. 

The Paper is organized as follows. Computational de¬ 
tails and the structural models are presented in Sec. [Ill 
In Section ED we discuss our results: structural effects 
(IIII AD . magnetoelectric coupling (IIIIBp . barrier metal¬ 
ization (inroii . and spin-dependent transport (HTTPl) . We 
close with a summary and an outlook (Sect. HVj) . 

II. STRUCTURE MODELS AND 
COMPUTATIONAL DETAILS 

The C0/PTO/C0 junctions are modeled by supercells 
(Fig. [T|) . In analogy to Co films on GaAs substrates^, we 
assume that the Co atoms form a face-centered tetrago¬ 
nal (fct) structure, which is in agreement with epitaxial 
growth of Co on PTO. 
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Because the barrier consists of complete PTO unit cells 
(UCs), in analogy to BaTi0 3 (Refs. [13 and |ll|) , we are 
concerned with two different interfaces with the Co elec¬ 
trodes on either side of the PTO stack: Co 2 /PbO (left 
interface) and Ti02/Co2 (right interface). In agreement 
with earlier studies of Ti 02 /Co 2 interfaces^, the most 
stable configuration has Co placed in line with the in¬ 
terfacial O ions. At a Co 2 /PbO interface, Co atoms are 
placed in line with Pb and O. 

The barrier thickness is five perovskite AB 0 3 UCs, 
the in-plane lattice constant is taken from experiment 
(PTO: a = 3.892 A, Ref. Eut ). At each side of the barrier 
five atomic layers of fct Co as well as a vacuum region 
of 20 A are attached. We define majority and minority 
spins with respect to the right Co electrode. 

To mimic a PZT barrier, we rely on a digital alloy 
model in which interfacial Ti is replaced by Zr. Thus, Zr 
forms a chemically complete Zr 02 monolayer at the right 
interface [Fig.QJb)]. Earlier studies^ showed that a simu¬ 
lation of the 25 % Zr composition in each Ti(Zr)02 layer 
by in-plane 2x2 supercells gives structural and magnetic 
properties of the Co/PZT interface which are similar to 
the digital alloy model used here. In the forthcoming, we 
refer to this barrier as PTO-ZO. 

The geometries of the above supercells have been op¬ 
timized using the VASP cod o 22 i 23 . Initially, the ionic dis¬ 
placements in the ferroelectric were taken as their theo¬ 
retical bulk values (0.45 A and 0.33 A for PbO and TiCU 
planes, respectively), in accordance with the ferroelec¬ 
tric polarization being along [001] (z axis). The fct Co 
leads have two atoms per layer with an interlayer distance 
of 1.47 A along [001] and an in-plane nearest-neighbor 
distance of 2.75 A. Within the barrier, the atomic po¬ 
sitions of the three central UCs were fixed to the bulk 
values but all atoms near the interfaces were allowed to 
move along the [001] direction. A T-centered 4x4x2 
Monkhorst-Pack mesh^i and an upper limit of 10 -2 eV/A 
for the ionic forces guarantee accurate structural relax¬ 
ations. We define P_^ for displacements A z > 0 and P<_ 
for Az < 0. 

For the transport calculations the above supercells 
have to be extended. To treat both the parallel (P, j"t) 
and antiparallel (AP, J/f) alignments of the magnetiza¬ 
tions within the Co electrodes, a double supercell was 
constructed by creating a mirror image of the original 
supercell with respect to the (001) plane and attaching 
it to the original junction [Fig. DJc)]. These double su¬ 
percells comprise two PTO barriers whose ferroelectric 
polarizations are oppositely oriented to each other. Co 
layers were added to ensure the correct layer alternation. 
The electronic states in the interior of the Co stacks have 
to be close to those of Co bulk since later, with respect 
to an appropriate treatment of electronic transport, two 
semi-infinite leads will be attached to the scattering re¬ 
gion. 

For the calculation of the tunnel conductances G, one 
half of a double supercell is taken as scattering region. 
Using the Landauer-Buttiker formalism implemented in 


the QUANTUM ESPRESSO code2£r— , transmission func¬ 
tions T(k x , k y ) have been computed for C 0 /PTO/C 0 and 
C 0 /PTO-ZO/C 0 within the two-dimensional Brillouin 
zone (2BZ). The latter is a square with —it/ a < k x < tt / a 
and —it/ a < k y < tt/ a. An adaptive k mesh refinement^ 
reliably yields accurate tunnel conductances 

G = [ T{k x , k v ) dk 2 . (1) 

J 2BZ 

Conductances with an error less than 10 -8 e 2 /h are 
achieved after a few refinements. To resolve all impor¬ 
tant features of the transmission maps, about 48000 k 
points in the 2BZ prove sufficient. The smallest trian¬ 
gles in the adaptive integration have an area of about 
5 • 10- 6 (27r/a) 2 . 

Having well converged tunnel conductances G, the 
TMR is calculated by 

rdt _ f4t 

TMR(P) = • 100 %, P = P<_, P-> ( 2 ) 

Gp + Gp 

for both orientations of the ferroelectric polarization P. 
For each magnetic configuration M of the leads, the TER 
ratio is given by 

gM _ gM 

TER(M) = + qm ■ 100 M =tt, U (3) 

in which G^f (G^) corresponds to the conductance for 
the barrier polarization pointing towards the TiC> 2 /Co 2 
(Co 2 /PbO) interface. 

Spin-resolved densities of states (DOS) have been cal¬ 
culated for these systems using the quantum ESPRESSO 
code2£. We use a T-centered 15 x 15 x 1 mesh with 
a smearing of 0.02 Ryd. Energy cutoffs read 63Ryd for 
the wavefunctions and 504 Ryd for the charge density. 
All calculations were performed within the generalized- 
gradient approximation^ (GGA-PBE) to the exchange- 
correlation potential. 

III. DISCUSSION AND RESULTS 

A. Structural effects 

The ionic displacements in the PbO and Ti02 planes 
of relaxed PTO are homogeneous along [001] for the 
P<_ configurations, that is, they do not depend signif¬ 
icantly on the layer. For the opposite polarization P_^, 
these displacements are considerably suppressed near the 
Ti02/Co2-terminated interface^, which is related to the 
lower stability of this ferroelectric configuration. Zr im¬ 
purities introduced into the PTO barrier stabilize the lo¬ 
cal displacements across the entire barrier; the displace¬ 
ment at the ‘digital’ Zr 02 layer is sizably increased, in 
agreement with a previous study of Co/PZT interfaces^. 

The effective thickness d e ^ of the ferroelectric barrier is 
an important property as far as transport properties are 
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FIG. 1: (Color online) Structural models for Co/PTO/Co (a) and Co/PTO-ZO/Co (b) multiferroic tunnel junctions for the 
two ferroelectric polarizations P<_ and P_>. of the barriers. The chemical compositions of the interfaces are marked below each 
structure. The effective thickness d e ff of the barrier, i. e. the distance between the Co electrodes, is indicated in (a). In (c), 
both an original (left, used in the geometry optimization) and a doubled supercell (right, used in the transport calculations) 
are sketched. 


concerned. It may simply be defined as the average dis¬ 
tance between the terminating Co layers of each electrode 
[Fig. DJa)]. We find that d e ff is reduced for P_^ in com¬ 
parison to P<_: 1.1 A for PTO and 0.8 A for PTO-ZO. 
This observation is in agreement with recent transport 
measurements of Co/PZT/LSMO tunnel junctions^. It 
is explained by the bonding at the interfaces. Whereas 
the Co-0 and Ti(Zr)-Co bonds dominate at the B O2/C02 
side (P = Ti, Zr), the Pb-Co and Co-0 chemical bonds 
determine the geometry of the Co2/PbO interface. In 
both cases, the ferroelectric switching leads to either con¬ 
traction or expansion of these bonds, thereby changing 


B. Magnetoelectric coupling 

The geometrical changes that accompany the reversal 
of the ferroelectric polarization affect the magnetic struc¬ 
ture at the interfaces as well. For a Co/PTO-ZO/Co 
MFTJ we find sizable magnetic moments induced on the 
O ions on the left side and on the Zr cations on the right 
side in comparison with those of the Co leads (1.8/ip; 
Fig. [2]). This behavior is present also in Co/PTO/Co 
junctions, for which even larger magnetic moments are 
induced on the Ti cations at the Ti02/Co2 interface. 
These findings indicate hybridization of orbitals at the 
interfaces, in particular those involved in the Ti-Co and 
Co-0 bonds^. As a result, a strong magnetoelectric (ME) 
coupling is found for both interfaces. 

A similar mechanism for the ME coupling is expected 



FIG. 2: (Color online) Magnetoelectric coupling in a 
Co/PTO-ZO/Co tunnel junction. Top: difference of the 
layer-resolved magnetic moments Mi upon reversal of the fer¬ 
roelectric polarization in the barrier (P«_ and P_>): A Mi = 
Mi(P->) — Mi(P«_), i layer index. Bottom: spin-resolved 
charge densities for the two ferroelectric polarizations (ma¬ 
jority: red; minority: blue). 


for the Co2/PbO interface. However, in this case, the 
effect is attributed mostly to the Co-0 bonds. As the 
polarization is switched from PV- to P^, the Co-0 bonds 
are shortened by more than 30 % (from 2.65 A to 1.78 A). 
The Co-Pb bonds, on the other hand, are moderately 
stretched from 2.51 A to 2.65 A. The combined effect 
of these structural changes is a buckling of the interfacial 
Co layers. Accordingly, the magnetic moments of the two 
Co species differ: the magnetic moments of Co atoms in 
line with Pb cations are considerably decreased by about 
—0.3/iB for P<_; simultaneously, those of the neighbor- 
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ing Pb cations are increased by about 0.2 The Co 
magnetic moments close to the interface vary between 
1.5 and 1.9 //b, depending on layer and ferroelectric 
polarization. 

The above effect is due to an increased overlap be¬ 
tween the electronic states of these cations; confer the 
site-projected density of states in Fig. 01 The magnetic 
moment of Co in line with O ions is less affected. How¬ 
ever, O-p orbitals hybridize with Co-d orbitals in a broad 
energy range. Furthermore, the charge density becomes 
more delocalized on the O sites for P_^ (top in Fig. [3]). 
The larger orbital overlap of the aforementioned states 
enhances the induced moments on the O sites which are 
parallel to those of the neighboring Co atoms. Therefore, 
the ME coupling at the Co 2 /PbO interface relies on the 
induced magnetic moments in the PbO layer and on the 
change of the Co magnetic moments. This strong local¬ 
ization of the ME coupling to the interface is evident from 
the layer-resolved magnetization in the Co/PTO-ZO/Co 
heterostructure (Fig. [2]). 

The strength of the ME coupling is quantified by the 
magnetoelectric coupling constant a = AM/(AE C ). It is 
defined as the change AM of the magnetization in the 
interface area A that is induced by the coercive electric 
field E c (that is the minimum field strength needed to 
switch the ferroelectric polarization). From the data pre¬ 
sented in Fig. [2j we estimate a to 9.1 • 10 -10 Gcm 2 V -1 
for the Co 2 /PbO and to —5.7 • 10 -10 Gcm 2 V -1 for the 
Ti 02 /Co 2 interface, assuming the coercive field of bulk 
PTO (E c = 33kV/cm; the coercive field for the PTO- 
ZO barrier might be smaller). These numbers are in 
line with those reported for similar heterostructures^i 
(e.g. Co/PZT, Ref.S LSMO/PTO, Ref.Ql Fe/BaTi0 3 , 
Ref. HI). 

Summarizing at this point, a strong ME coupling of 
electronic origin is established for both interfaces in the 
undoped and in the Zr-doped Co/PTO/Co junctions. 
The two interfaces differ with respect to the sign of 
the magnetization change in response to the ferroelectric 
switching. We find no qualitative change upon Zr dop¬ 
ing of the PTO barrier with respect to the undoped bar¬ 
rier. The ME coupling is traced back to the polarization- 
dependent hybridization of orbitals that are involved in 
the Co-O, Pb-Co, and Ti(Zr)-0 bonds. Because the de¬ 
tailed electronic structure at the interfaces determines 
significantly the spin-dependent transport, we expect a 
considerable effect on the electron transmission across the 
junctions. 


C. Barrier metalization 

Before discussing spin-dependent transport, we ad¬ 
dress whether the barriers are insulating (tunneling 
regime) or conducting (metallic regime). For this pur¬ 
pose, we investigate the site-resolved density of states at 
the Fermi level; a layer is considered insulating if this 
quantity is negligibly small. 




E - E f (eV) 


FIG. 3: (Color online) Effect of the ferroelectric polarization 
P in a Co/PZT/Co tunnel junction on the electronic struc¬ 
ture at the interfaces for (top) and P<_ (bottom) configu¬ 
rations. Spin-resolved densities of states are depicted for Co 
atoms in line with Pb [Co (Pb)] and in line with O [Co (O)], 
as well as for O and Pb cations at the Co 2 /PbO interface. 


For the P<_ configuration, the PTO barrier is locally 
metallic within one UC at the Co 2 /PbO interface and in 
the Ti02 plane at the right interface; the inner 3.5 UCs 
are insulating [Fig. 0(a)]. This partial metalization is due 
to charge transfer from Co atoms into the ferroelectric, 
which happens at both interfaces. For P_^, the effect is 
even enhanced, so that the entire barrier becomes metal¬ 
lic [Fig. 0(b)] and one expects very large conductances 
(as compared to those for the tunneling regime). 

The metalization affects mostly the Ti02 layers but 
leaves the PbO layers insulating. An exception is the 
terminating PbO layer at the left interface, in which siz¬ 
able magnetic moments are induced on the O sites. Such 
an alternating sequence of conducting Ti 02 and insulat¬ 
ing PbO layers might open new possibilities for electron 
transport. 

In contrast to the cases discussed above, the Zr-doped 
barriers are insulating irrespectively of the ferroelectric 
polarization [Figs. 0(c) and (d)] and, thus, we are con¬ 
cerned with the tunneling regime. More precisely, metal- 
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FIG. 4: (Color online) Electronic structure of Co/PTO/Co [upper two panels, (a) and (b)] and Co/PTO-ZO/Co [lower two 
panels, (c) and (d)] tunnel junctions for parallel alignment of the lead magnetizations (black vertical arrows, tt) an d for the 
two polarizations [(a) and (c)] and [(b) and (d)] in the ferroelectric barrier (as indicated in each panel by horizontal 
arrows). Spin-resolved densities of states are shown for the series of unit cells across the barriers; majority spin red, minority 
spin blue. 
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TABLE I: Spin-resolved conductances G of Co/PTO/Co mul- 
tiferroic tunnel junctions (in units of 10 _6 e 2 /h) for the par¬ 
allel (tt) and the antiparallel (4%) magnetic configurations of 
the Co electrodes as well as for the two ferroelectric polar¬ 
ization of the barriers (P_^ and P<_). Minority (min.) and 
majority (maj.) spin channels are defined with respect to the 
right Co electrode. The respective tunnel magnetoresistances 
(TMR) and tunnel electroresistances (TER) are given in ad¬ 
dition. 


PTO 

P^ 

P-+ 



Cmin Cmaj 

Cmin Cmaj 

TER 

tt 

1.40 2.03 

159 11.3 

96.1% 

tt 

0.82 0.80 

15.7 117.8 

97.6% 

TMR 

36% 

12% 



ization is restricted to a single UC at each interface, simi¬ 
larly to Co/PTO/Co for P«_ [Fig. Ufa)]. It is conceivable 
that this fundamental difference of the PTO- and PTO- 
ZO-based junctions is correlated with the aforementioned 
stabilization of the ionic displacements in the barrier by 
Zr impurities. 

D. Spin-dependent transport 

We now discuss spin-dependent transport for the two 
types of tunnel junctions in four different configurations. 

1. Co/PTO/Co 

For the Co/PTO/Co tunnel junction, the transmit- 
tances T(k x ,k y ) reveal a drastic change of the transmis¬ 
sion (upper part of Fig.[5j): for P<_ the ferroelectric PTO 
barrier is insulating, whereas for P_^ it is conducting, as 
shown in Fig.SJ In the latter case, large parts of the 2BZ 
show very high transmittances (red regions). 

The transmittance maps for P<_ (upper left part of 
Fig. [5]) manifest that the parallel magnetic configuration 
provides larger conductances than the antiparallel con¬ 
figuration, as is also evident from the numbers given in 
Table ID Thus, we are concerned with a normal rather 
than with an inverted TMR. 

For the parallel configuration of the electrodes (tt)> 
transport is dominated by free-electron-like majority 
electrons. The largest transmission is observed near the 
center of the 2BZ. The minority channel also contributes 
to the transmission but shows a complicated behavior, 
with the largest transmission at the X [fc(X) = (^/a, 0) 
and equivalent] and at the M [fc(M) = (^/a^/a) and 
equivalent] points of the 2BZ. The antiparallel config¬ 
uration (J%) is characterized by reduced transmittances. 
Here, both spin channels contribute almost equally to the 
conductance. 

Turning to the P_^ configuration (upper right part of 
Fig. [5j), the transmission maps change dramatically; we 


TABLE II: As Table [I] but for Co/PTO-ZO/Co multiferroic 
tunnel junctions. 


PTO-ZO 

p «- 

P ■-> 



G m in C ma j 

G min C ma j 

TER 

tt 

1.43 3.55 

16.2 0.79 

55% 

tt 

0.70 0.82 

0.80 8.07 

71% 

TMR 

53% 

31% 



recall that here the barrier becomes conducting, lead¬ 
ing to sizable areas in the 2BZ with high transmittance. 
There are also large contributions that arise from res¬ 
onant tunneling. As a result, the conductances are in¬ 
creased by two orders of magnitude as compared to 
the tunneling regime, yielding giant TER ratios (almost 
100%, Table |J). 

It turns out that also for the metallic PTO barrier the 
TMR is not inverted, in contrast to experiment. A TMR 
ratio of 12 % is significantly lower than for P^_ (36 %). 

We conclude that—although being illustrative with re¬ 
spect to effects of insulating and conducting barriers— 
the PTO barrier cannot model reliably the experimental 
observations. 


2. Co/PTO-ZO/Co 

The effect of Zr doping on the electronic trans¬ 
port is analyzed by comparing the transmittances of 
Co/PTO/Co with those of the Co/PTO-ZO/Co digital 
alloy. 

For P^_ (bottom left part of Fig. [5]), the transmission 
for the parallel configuration is enhanced in the majority 
spin channel for the Zr-doped system, whereas it is mildly 
affected in the minority spin channel (cf. Table HU). This 
observation is at variance with an admittedly very simple 
explanation by the effective thickness of the ferroelectric 
barrier; the latter increases upon Zr doping, which would 
usually lead to lower transmittances because of the expo¬ 
nential decay of electronic states within the barrier. One 
could argue that the larger extent of the Zr-d orbitals 
compared to the Ti-d orbitals might increase the trans¬ 
mission probability across the Zr 02 /Co 2 interface due to 
the stronger hybridization with the orbitals at neighbor¬ 
ing atomic sites. 

A significant change of the conductances upon Zr dop¬ 
ing is not observed for the antiparallel lead magnetiza¬ 
tions (Table ITT]) . In this case, the Zr impurities result in a 
redistribution of the transmittances but keep the general 
shape of the T maps. 

The transmittance maps for P_^ display an eye¬ 
catching difference to those for P<_. Nevertheless, the 
sign of the TMR is preserved and its absolute value is 
moderately reduced from 53 % to 31 %. We recall that 
for P_^ we are concerned with tunneling, in contrast to 
the Co/PTO/Co junction. 
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FIG. 5: (Color online) Spin-dependent transport in Co/PTO/Co (top row) and Co/PTO-ZO/Co (bottom row) tunnel junctions. 
In each panel, transmission maps T(k x , k y ) are shown as a color scale for the parallel (Tt) and the antiparallel (]/[) configuration 
of magnetic Co electrodes as well as for the two ferroelectric polarizations of the barriers (indicated by horizontal arrows). Color 
scales in units of 10“ 5 e 2 /h. 


At the Fermi level, the electronic properties of Co 
are mainly determined by the majority-spin s and the 
minority-spin d states. Although the minority states 
dominate in bulk Co at the Fermi level, our ab initio cal¬ 
culations for P<_ tunnel junctions suggest that these elec¬ 
trons are transmitted less than the majority electrons. 

Concerning P^, the electron (hole) charge transfer 
across the interface plays a crucial role. We found that 
a large (in absolute value) negative spin polarization is 
induced on the interfacial Ti/Zr-d orbitals, which origi¬ 
nates from the hybridization with the neighboring d or¬ 
bitals of Co (cf. Figs. 11 and 12 in Ref. |9|). As a result, 
the by far largest conductances show up in the minority 
channel for tt and in the majority channel for tt (Ta¬ 
ble HU). 


The TER is explained by the distribution of free charge 
carriers in this heterostructure. According to the layer- 
resolved DOS (Fig. |4]), an appropriately oriented ferro¬ 
electric polarization of the barrier induces charge transfer 
at the Ti 02 /Co 2 interface which can lead to enhanced 
propagation of electronic states. This effect would re¬ 
duce the decay length of electrons within the barrier for 
P_^ as compared to P<_. Since the tunneling probability 
increases exponentially with decrease of the barrier thick¬ 
ness, the above-mentioned partial metalization leads to 
larger total conductances for P_^. This scenario is con¬ 
firmed by our ab initio results (Fig. [5] and Table HU) . 
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IV. SUMMARY 

The present ab initio study clarifies the origin of the 
tunnel magnetoresistance (TMR) and the tunnel elec¬ 
troresistance (TER) in multiferroic Co/PbTiOs/Co. The 
results may be ‘transferred’ to similar heterostructures, 
i. e. those that combine ferromagnetic electrodes (e. g. 
Fe, Co, and Ni) and well-established ferroelectrics (e. g. 
BaTiOs). The strong magnetoelectric coupling at the 
Co/PbTi 03 interfaces is explained by spin-dependent hy¬ 
bridization of orbitals; it affects considerably the local 
magnetic moments and the spin-resolved charge density 
in a narrow region about the interface. As a result, we 
find four distinctly different conductances, that is a four- 
state memory device. 

The charge transfer from Co into the ferroelectric bar¬ 
rier results in a metalization of the latter. This trans¬ 
fer is particularly strong for the ferroelectric polarization 
pointing towards the Ti02/Co2 interface. Such a com¬ 
plete metalization may be attributed to the local-density 
approximation in the underlying density-functional cal¬ 
culations which show up as a too small fundamental band 
gap (1.7 eV in the bulk but about 1.5 eV within the bar¬ 
rier; experiment: 3.4 eV, from Ref. [33] and references 
therein): if the calculated band gap would approach the 


experimental gapi^, the insulating region of the PbTiOs 
barrier would be almost as wide as its nominal chemi¬ 
cal thickness. This issue needs to be clarified in a future 
investigation. 

The inverted TMR observed in experiments on 
Co/PbZro^Tio.sOs/LaSrMnOs tunnel junctions^ is not 
reproduced by calculations for the digitally Zr-doped 
PbTiOs barrier sandwiched between Co electrodes. This 
finding suggests that La(Sr)Mn 03 plays a crucial role 
for the TMR inversion. Within this respect, the almost 
half-metallic MnC^/PbO interface, instead of Co 2 /PbO 
considered in this study, may significantly modify the se¬ 
lection of states for spin-dependent transmission (‘sym¬ 
metry filtering ,>34-36. 

Regarding the PbZro. 2 Tio. 8 O 3 bar¬ 
rier, randomly distributed Zr dopants may influence the 
tunneling as well, since Zr enhances locally the ferroelec¬ 
tric displacements in PbTiOs and modifies the bonding 
at the interfaces^. These issues will treated in a future 
theoretical study. 
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